ON THE ASYMPTOTICS OF A TOEPLITZ DETERMINANT WITH 

SINGULARITIES 



P. DEIFT, A. ITS, AND I. KRASOVSKY 

Abstract. We provide an alternative proof of the classical single-term asymptotics for Toeplitz 
f*^ , determinants whose symbols possess Fisher-Hartwig singularities. We also relax the smoothness 

■ conditions on the regular part of the symbols and obtain an estimate for the error term in the 

asymptotics. Our proof is based on the Riemann-Hilbert analysis of the related systems of orthogo- 
nal polynomials and on differential identities for Toeplitz determinants. The result discussed in this 
paper is crucial for the proof of the asymptotics in the general case of Fisher-Hartwig singularities 
and extensions to Hankel and Toeplitz+Hankel determinants in [15]. 
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Let f(z) be a complex- valued function integrable over the unit circle. Denote its Fourier coeffi- 
cients 

1 f 2n 

fj = — f(e i6 )e-^de, j = 0,±l,±2,... 
27T Jo 

We are interested in the n-dimensional Toeplitz determinant with symbol f(z), 
(1-1) D n {f{z)) = det(/ J -_ fc )^ 0) n > 1, 

where f{e ie ) has a fixed number of Fisher-Hartwig singularities [20, 27], i.e., / has the following 
form on the unit circle: 
*sO '. m 
& i (1.2) f(z) = e v ^z^o? 3 Y[\z- Zj \ 2 ^g ZjA (z)z/\ z = e w , 9 € [0,2tt), 

j=0 

for some m = 0, 1, . . . , where 
(1.3) Zj = e i9 i, j = 0,...,m, = 9 <9 1 <--- < 6 m < 2tt; 

I e Mj < argz < 2-k 
(1.5) SRqj- > -1/2, 0j€C, j = 0,...,m, 

and y(e iS ) is a sufficiently smooth function on the unit circle (see below). The condition on the 
ctj 's insures integrability. Note that a single Fisher-Hartwig singularity at Zj consists of a root-type 
singularity 



(1.6) \z - Zj 



\2cLj 



2 sin ■ 



and a jump e J7r/3 — > e~™$ . We assume that Zj, j = 1, . . . ,m, are genuine singular points, i.e., either 
ay ^ or /3j ^ 0. However, we always include Zq = 1 explicitly in (1.2), even when ao = /9o = 0: 
this convention was adopted in [15] in order to facilitate the application of our Toeplitz methods 
to Hankel determinants. Note that gp {z) = e _l7r/3 °. Observe that for each j ^ 0, z^g^.{z) is 
continuous at z = 1, and so for each j each "beta" singularity produces a jump only at the point 

l 
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Zj. The factors z- 3 are singled out to simplify comparisons with the existing literature. Indeed, 
(1.2) with the notation 6(0) = e v ^ is exactly the symbol considered in [20, 3, 4, 8, 9, 10, 18, 19, 29]. 
However, we write the symbol in a form with z^= oPj factored out. The representation (1.2) is 
more natural for our analysis. 

On the unit circle, V(z) is represented by its Fourier expansion: 

00 1 /-27T 

(1.7) V( Z )=J2 V ^^ Vk = ^J V ( ei6 ) e ~ ki9de - 

k=—oo ® 

The canonical Wiener-Hopf factorization of e v ^ is given by 

(1.8) e v W =b+(z)e Vo b-(z), b+(z) = e^=i v k *\ = e T. k L x v k z\ 
Define a seminorm 

(1.9) III^HI = max 

J,k 

where the indices j, k = are omitted if z = 1 is not a singular point, i.e. if ao = Po = 0. If m = 0, 
set = 0. 

In this paper we consider the asymptotics of D n (f), n — > oo, in the case \\\/3\\\ < 1. The case 
\\\/3\\\ > 1 was addressed in [15]. The asymptotic behavior of D n (f) has been studied by many 
authors (see [15, 19] for a review). An expansion of D n (f) for the case V € C°°, < 1, was 

obtained by Ehrhardt in [19]. The aim of this paper is to provide an alternative proof of this 
result based on differential identities for D n (f) and a Riemann-Hilbert-Problem analysis of the 
corresponding system of orthogonal polynomials. We also obtain estimates for the error term and 
extend the validity of the result to less smooth V. The analysis of D n {f) for |||/3||| < 1 plays a 
crucial role in the analysis of D n (f) for > 1 (see [15]). 

We prove 

Theorem 1.1. Let f{e ie ) be defined in (1.2), |||/3||| < 1, Ray > -1/2, ay ± pj + -1,-2,... for 
j, k = 0,1, . . . ,m, and let V(z) satisfy the condition (1.11), (1.12) below. Then as n — > oo, 



(1.10) D n (f) = exp 



nV + J2 kVkV - 



k=l J j=0 

aj/3 k -a k /3j 



x „£r=o(°?-^) J] \ Zj - Zk \2(^ k -a jak ) 



0<j<k<rn 

m 



T-r Gq + aj + ^Ml + aj-^) 

n g(Tt^) (1 + o(1)) ' 



where G{x) is Barnes' G-function. The double product over j < k is set to 1 ifm = 0. 

Remark 1.2. As indicated above, this result was first obtained by Ehrhardt in the case V € C°°. 
We prove the theorem for V(z) satisfying the smoothness condition 

oo 

(i.ii) Yl i*mi<°°, 

k=— oo 

where 

ni2 , s , i + E.-Lo^^ + TO) 2 ] 

[ ' i- HI/8111 
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Remark 1.3. In the case of a single singularity, i.e., when m = 1 and ocq = /?o = 0, or when m = 0, 
the seminorm |||/3||| =0, and the theorem implies that the asymptotic form (1.10) holds for all 

(1.13) Ka m >-i, /3 m €C, a m ±/3 m / -1,-2,... 

In fact, if there is only one singularity, say m = 0, and V = 0, an explicit formula was found by 
Bottcher and Silbermann [8] for D n (f) for any n in terms of the G-functions: 

G(l + a + f3 )G(l + a - go) G(n + l)G(n + 1 + 2a ) 

^ 114 j ™ l j G(l + 2a ) G(n + l + a + /5o)G , (n + l + a -/5o)' 

^cto > °o ± A) / -1, -2, . . . , n > 1, 

and (1.10) can then be read off from the known asymptotics of the G-function (see, e.g., [2]). 

Remark 1.4. Assume that the function V(z) is sufficiently smooth, i.e. such that s in (1.11) is, in 
addition to satisfying (1.12), sufficiently large in comparison with |||/3|||. Then we show that the 
error term o(l) = ©(n'"^"' -1 ) in (1.10). In particular, the error term o(l) = 0(ri' 1 1' 9 ! 1 if V(z) 
is analytic in a neighborhood of the unit circle. Moreover, for analytic V(z), our methods would 
allow us to calculate, in principle, the full asymptotic expansion rather than just the leading term 
presented in (1.10). Various regularity properties of the expansion (uniformity, differentiability) in 
compact sets of parameters satisfying 9ftay > —1/2, |||/3[|| < 1, ctj ± (3j ^ —1, —2, ... ,6j ^ 8^, are 
easy to deduce from our analysis. 

Remark 1.5. Since G(—k) = 0, k = 0, 1, . . . , the formula (1.10) no longer represents the leading 
asymptotics if ctj + /3j or ctj — j3j is a negative integer for some j. Although our method applies, we 
do not address these cases in this paper. It is simply a matter of going deeper in the asymptotic 
expansion for D n (f). It can happen that D n (f) vanishes to all orders (cf. discussion of the Ising 
model at temperatures above the critical temperature in [16]), and e~ nVo D n (f) is exponentially 
decreasing. 

We prove Theorem 1.1 in the following way. We begin by deriving differential identities for the 
logarithm of D n (f) in Section 3, in the spirit of [17, 25, 26, 22, 13, 14], utilizing the polynomials 
orthogonal with respect to the weight f(z) on the unit circle. Then, assuming that V{z) is analytic 
in a neighborhood of the unit circle, we analyze in Section 4 the asymptotics of these polynomials 
using Riemann-Hilbert/steepest-descent methods as in [15]. This gives in turn the asymptotics of 
the differential identities from which the formula (1.10) follows in the V = case by integration 
w.r.t. ctj, f3j in Section 5.1. However, the error term that results is of order n 2 !!^!!" 1 Inn (see 
(5.35)), which is asymptotically small only for |||/3||| < 1/2, rather than in the full range |||/3||| < 1. 
To prove (1.10) for all \\\/3\\\ < 1, we need a finer analysis of cancellations in the Riemann-Hilbert 
problem as n — > oo. We carry this out in Section 5.2 and reduce the leading order terms in D n (f) 

to a telescopic form (see (5.65)), which leads to a uniform bound on D n (f)n^j= ^j~ a ^ for large 
n which is valid for all |||/3||| < 1 and away from the points ctj ± f3j = —1, —2, .... We then apply 
Vitali's theorem together with the previous result for \\\/3\\\ < 1/2. This proves Theorem 1.1 in the 
V = case as desired, with the error term of order n'"^'" _1 . In Sections 5.3, 5.4, we then extend 
the result to the case of analytic V ^ by applying another differential identity from Section 3. 

Ehrhardt proves (1.10) using a "localization" or "separation" technique, introduced by Basor in 
[4], in which the effect of adding in Fisher-Hartwig singularities one at a time, is controlled. One 
may also view our approach as a "separation" technique, but in contrast to [4, 19], we add in the 
Fisher-Hartwig singularities, as well as the regular term e v ( z \ in a continuous fashion. 
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Remark 1.6. An alternative approach to proving Theorem 1.1 in the ^-analytic case is to apply, ab 
initio, the finer analysis of Section 5.2 to the orthogonal polynomials which appear in the differential 
identities. This approach is more direct, but is considerably more involved technically. The analysis 
is resolved, as above, by reduction of the problem to an appropriate telescopic form. 

Finally in Section 5.5, we extend our result to the case when V(z) is not analytic and only satisfies 
the smoothness condition (1.11), (1.12). We approximate such V{z) by trigonometric polynomials 
V^ a \z) = Yl P k=- p Vk zk with an appropriate p = p{n) and modify the Riemann-Hilbert analysis 
accordingly. This produces asymptotics of a Toeplitz determinant in whose symbol f^ n \z) the 
function V(z) is replaced by V^ n \z). We then use the Heine representation of Toeplitz determinants 
by multiple integrals to show that D n (f^) approximates D n (f) as n — > oo, sufficiently strongly 
to conclude Theorem 1.1 in the general case. 

2. Riemann-Hilbert problem 

In this section we formulate a Riemann-Hilbert problem (RHP) for the polynomials orthogonal 
on the unit circle (which, oriented in the positive direction, we denote C) with weight f{z) given 
by (1.2). We use this RHP in Section 4 to find the asymptotics of the polynomials in the case of 
analytic V(z). Suppose that all Dk{f) ^ 0, k = ko, ko + 1 . . . , for some sufficiently large ko (see 
discussion below). Then the polynomials 4>k(z) = XkZ k + ■ ■ ■ , 4>k(z) = XkZ k + ■■• of degree k, 
k = ko, ko + 1, . . . , satisfying the orthogonality conditions 



(2.1) 



1 

27 



2tt 



<j> k (z)z-*f(z)dO 



Xk^jk, 



1 

27 



2tt 



<t> k {z~ L )zP f{z)dB 



X k ^jk, 



0,1,..., k, 



exist and are given by the following expressions: 

/oo 

(2.2) 4> k (z) = 



\/Dk~Dk] 



/fe-10 

1 



foi 
hi 

fk-ll 

Z 



fok 
flk 



fk-lk 
^k 



(2.3) 



where 



\f~DkP, 



fc+i 



/oo /oi 

/lO /ll 



/ofe-i 1 

/lfc-1 z~ l 



fkO fkl ••• fkk-1 Z k 
1 t' ZlT 

/st= 2W f^ z ~ {S ~ t)dd > s,t = 0,l,...,k. 



We obviously have 
(2.4) 



Xk 



D 



k+i 



(2.5) 



Consider the following 2x2 matrix valued function Y ( ~ k \z) = Y(z), k > ko'- 

x k x 4>k{z) Xk 1 Ic 



y(fe)( z ) = 



-ij. f ^ , -i r MO fij}S \ 

IC €-z 2ni£ k \ 

-Xk-^-^k-iiz' 1 ) -Xk-i Jc ?fc - l(rl)/m I ' 
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It is easy to verify that Y(z) solves the following Riemann-Hilbert problem: 

(a) Y{z) is analytic for z € C \ C. 

(b) Let z £ C\yj 1 J l = QZj. Y has continuous boundary values Y + (z) as z approaches the unit circle 
from the inside, and Y"_(z), from the outside, related by the jump condition 

(2.6) Y + {z) = Y_{z)(\ zeC\Uf =oZj . 



(2.7) ) I I + ()(}-))(-* ), ;w: . 



(c) has the following asymptotic behavior at infinity: 

z k ( 

(d) As z — >• Zj, j = 0, 1, . . . ,m, z G C \ C, 

f2 8) Y(z)-(° {1) 0{l) + 0{\z-z 3 \^)\ 

(2-8) r ^"^0(l) O(l) + O(|z-z,f a 0y ' lta ^°' 



and 



(2.9) 



y(z) _fO(l) 0(ln|.-^|)\ 



(Here and below O(a) stands for 0(|a[).) If ao = A) = 0, F(z) is bounded at z = 1. 

A general fact that orthogonal polynomials can be so represented as a solution of a Riemann- 
Hilbert problem was noticed in [21] for polynomials on the line and extended to polynomials on the 
circle in [1]. This fact is important because it turns out that the RHP can be efficiently analyzed for 
large k by a steepest-descent type method found in [11] and developed further in many subsequent 
works. Thus, we first find the solution to the problem (a)-(d) for large k (applying this method) 
and then interpret it as the asymptotics of the orthogonal polynomials by (2.5). 

Recall the Heine representation for a Toeplitz determinant: 

(2.10) d m=jmJ -J n i^- e ^i 2 n/(^)#,- 

If f(z) is positive on the unit circle, it follows from (2.10) that Df t (f) > 0, k = 1, 2, . . . , i.e. ko = 1. 
In the general case, let A be a compact subset in the subset |||/3||| < 1, ay ± /3j / —1,-2,... 
of the parameter space V = {(ao, /3o, . . . , a m , (3 m ) : ctj,/3j € C, Kay > —1/2}. We will show in 
Section 4 that the RHP (a)-(d) is solvable in A, in particular \k are finite and nonzero, for all 
sufficiently large k (k > fco(A)). Let fifc be the set of parameters in V such that -Dfc(/) = for 
some k = 1, 2, . . . , ko — 1. We will then have Dk(f) ^ 0, k = 1, 2, . . . for all points in A \ Ofc . 
Note that D n (f) depends analytically on ay, (3j in V: this is true, in particular, on (the interior 

of) A\n h) . 

The solution to the RHP (a)-(d) is unique. Note first that detY(z) = 1. Indeed, from the 
conditions on Y(z), det Y(z) is analytic across the unit circle, has all singularities removable, and 
tends to 1 as z —> oo. It is then identically 1 by Liouville's theorem. Now if there is another 
solution Y(z), we easily obtain by Liouville's theorem that Y(z)Y(z)~ l = 1. 

3. Differential identities 

In this section we derive expressions for the derivative (d/d'y) In D n (f(z)), where either 7 = ctj 
or 7 = Pj, j = 0, 1, . . . , m, in terms of the matrix elements of (2.5). These will be exact differential 
identities valid for all n = 1, 2, . . . (see Proposition 3.1 below), provided all the D n (f) ^ 0. We will 
use these expressions in Section 5.1 to obtain the asymptotics (1.10) in the case V = 0, |||/?||| < 1/2 
(improved to |||/3||| < 1 in Section 5.2). Furthermore, in this section we will derive a differential 



6 P. DEIFT, A. ITS, AND I. KRASOVSKY 

identity (see Proposition 3.3 below) which will enable us in Sections 5.3, 5.4 to extend the results 
to analytic V ^ (improved to sufficiently smooth V in Section 5.5). 

Set Dq = 1, (f>o(z) = 4>o(z) = 1, and suppose that D n (f) ^ for all n = 1, 2, . . . . Then the 
orthogonal polynomials (2.2), (2.3) exist and are analytic in the Oj's, /3j's for all k = 1,2,.... 
Moreover, (2.4) implies that 

n— 1 

(3-1) D n (f(z)) = Ux7 2 



j=0 



Note that by orthogonality, 

f27T 



i r 



— y y"^"^ ~*~ P or y nonna l °f degree j — lj <f>j(z~ 1 )f(z)d6 = — 

Similarly, 

(3.3) ±[*M'&i£±to)*- 10 » 



2vr 7 <9 7 Xj dl 

Therefore, using equation (3.1), we obtain 

n r> U— 1 n— 1 ^Xj ^ „27T Q / n — 1 \ 

(3.4) ^lnD n (f(z)) = ^Uxj 2 = " 2 E^ = I ^ ( E^C*^- 1 )) 

Using here the Christoffel-Darboux identity (see, e.g., Lemma 2.3 of [15]) 

n-i ^ ^ / ^ d ^ \ 

^ 0Jfc(2 _1 )0fc(z) = -TKpniz^niz' 1 ) + Z ( (f) n (z' l ) — (f) n (z) - <p n (z) — (j> n (z' 1 ) J , 
fc=0 V Z 2 / 

and then orthogonality, we can write 
,3.5) _ lnD „,/W) = 2n^ + -y o 



5 1 n fJV U O "^f , 1 /" 27r 9 ( A ( \ d ^n{z X ) 2 ( -U d( t>n(z) \ 
' - •'■ ■ 1 ' ' U ' 1 — ' — ' 9n{z) ~. (f)n{z ) Zf(z)d9. 



07 \ dz dz 

Writing out the derivative w.r.t. 7 in the integral and using orthogonality, we obtain: 

(3.6) ^\nD n (f(z)) = h-I 2 , 

where 

(3 - 7) /l = ^y — /2 = ^U — — f{z)de - 

It turns out that the particular structure of Fisher-Hartwig singularities allows us to reduce (3.6) 
to a local formula, i.e. to replace the integrals by the polynomials (and their Cauchy transforms) 
evaluated only at several points (cf. [22, 26]). 

Let us encircle each of the points Zj by a sufficiently small disc, 

(3.8) U Zj ={z: \z-Zj\ <£}. 
Denote 

(3.9) c £ = u™ o (u Zj n c) . 
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We now integrate I\ by parts. First assume that V(z) = 0. Then, using the expression 
df(z) 



de 



3=0 



= E a J cot — o- 1 + ) /(*) 



E 

v3=0 



2 + Z 3 
' Z — Zj 



we obtain: 



(3-10) h = -Xn^\n 




lim 



9<t>n(z) - 
C\C e #7 



^) E 



Z + Zo 



O ; 



,3=0 



Z — Zj 



f(z)de 



iE^^^'H ' • )- /u,< ,: » 



3=0 ' 



where the integration is over C\C £ in the positive direction around the unit circle. 
Note that by adding and subtracting d ^ >r ^ 3 - , and by using orthogonality we can write 



(3.11) 



/ 



d(j>n(z) 
C\C e #7 

d<t>n{z) 



Mz~ l ) 



Z + Zj 



Z — Zj 



f(z)de 



z — z 



d7 Jc\Q 



2z 



Mz' 1 ) 3 —f(z)d0 + O(e 



Z — Zj 



Obviously, the fraction in the first integral on the r.h.s. is a polynomial in z of degree n — 1 with 
leading coefficient dx-n/dj. Therefore, the integral equals 27r^ L /xn up to 0(e 2?Raj+1 ). The second 
integral can be written in terms of the element Y22 of (2.5) for z — > Zj. Let us estimate therefore 
the following expression for ctj / 0: 



(3.12) 



/ Ms' 1 ) 
Jc s 



2zjf(s) ds 



lim 

z is e->o 



z — > z 



■31 



\z\ > 1. 



This difference tends to zero as z — > Zj, for -Ray > 0. When Ray < 0, it is a growing function as 
z — s> Zj, and when Ray = 0, Say 7^ 0, an oscillating one. The analysis is similar to that of Section 
3 in [26]. For future use, we now fix an analytical continuation of the absolute value, namely, write 
for z on the unit circle, 



(3.13) 



3> V Z 3 ) {zz . e U j)aj /2> 



\Z-Zj\~> 



where ij is found from the condition that the argument of the above function is zero on the unit 
circle. Let us fix the cut of (z — Zj) a] going along the line 6 = 9j from Zj to infinity. Fix the branch 
by the condition that on the line going from Zj to the right parallel to the real axis, arg(z — Zj) = 2ir. 
For z q j/ 2 in the denominator, < argz < 2ir. If zq = 1 let < arg(z — 1) < 2ir. (This choice 
will enable us to use the standard asymptotics for a confluent hypergeometric function in the RH 
analysis in Section 4 below.) Then, a simple consideration of triangles shows that 



(3.14) 



3tt, < 9 < 0j 
7T, 9j <9 < 2tt 
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Thus (3.13) is continued analytically to neighborhoods of the arcs < 9 < 9j, and 9j < 9 < 2ir. 
We now analyze (3.12) in the same way that equation (26) was analyzed in [26]. For this analysis, 
however, we will need two other choices of the function (z — Zj) 2a ^: one choice with the cut going 
a short distance clockwise along the unit circle C from Zj, and another, with the cut going a 
short distance anticlockwise along C from zj. Let Cj and dj be some points on C between zj 
and the neighboring singularity in the clockwise and anticlockwise directions, respectively. In a 



neighborhood of Zj, let g(s) be defined by the formula <j) n {s 1 )f(s)/(is) = \s 
obtain as in [26] for the part of (3.12) on the arc (cj,dj): 



1 2a. 



g(s). We then 



(3.15) 

rd 



' \s-zA 2a ' 



q(s)ds — lim 
s- z w e->0 



/ • • ./ 



Zj \ 2a i 



e 2a 3 



g{s)ds - —(gizje t£ ) - g(z je t£ )) 



ir(z — Zj) J 
hm — - 

e^o \zzj) a i sin(27raj) 



{e 2nia i- ialR g{zje- i£ )-e 
+ aj 1 0(z-z j ), 



2a 



2Trictj—ia£L 



z — > Zj, \z\ > 1, 

for a.j ^ 0, 1/2, 1, 3/2, . . . (for a.j = 1/2, 1, 3/2, . . . one obtains terms involving [z — Zj) k ln(z — z 
vanishing as z — >■ Zj). Here the constants £r, £l depend on the choice of a branch for (z — Zj 
(whose cut is, recall, along the circle) and their values will not be important below. 
Introduce a "regularized" version of the integral in a neighborhood of Zf. 



■3) 

2a 



(3.16) 



d^) \ s _ z .\2 aj 

s — z 



g(s)ds 



— lim 



d ' Is - zA 2a i 



s — z 
tt(z — Zj) 2a l 



g(s)ds 



e^o [zzj) a i sin(27raj) 



2iriaj —iaiu 



ie\ —2-wiaj—ialL 



for z in a complex neighborhood of Zj and —1/2 < Kay < 0, a 3 ^ 0. If 5?Qj > 0, we set the 
"regularized" integral equal to the integral itself. 

Denote by Y the matrix (2.5), in which the integrals of the second column are replaced by their 
"regularized" values in a neighborhood of each zj. 

Then, collecting our observations together, we can write (3.10) in the form 



(3.17) h = - X „ 



n 



3=0 



3=0 



2a i z j d 



'-(.n+l). 



a,- 



where 
(3.18) 

A similar analysis yields for I2: 

-idx 



Af(z j ) = ton{f(z j e-")-f{z j e»)). 



(3.19) I 2 = Xn 



n 



+ E(°y-&) + E 



3=0 



3=0 



1 



3d-y 



x -^ +1 \z 3 ))Y^(z 3 ), aj^O 



(n), 



k 2-ni dy 



-n(Zj)Af( Zj ), 



a.; 



0. 



Substituting these results into (3.6) we obtain 
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Proposition 3.1. Let V(z) = 0. Let 7 = oik or 7 = (3k, k = 0, 1, . . . , m, and D n (f(z)) ^ for all 
n. Then for any n = 1, 2, . . . , 



(3.20) il-D,(/(z)) = - 2x ;'|! 



j"=o 



ctj = 0, 



where Af(zj) is defined in (3.18). 

In Section 5.1 we substitute the asymptotics for Y (found in Section 4) in (3.20) and, by in- 
tegrating, obtain part of Theorem 1.1 for f(z) with V(z) = 0, \\\(3\\\ < 1/2. Further analysis of 
Section 5.2 extends the result to \\\6\\\ < 1. 



Remark 3.2. The differential identities (3.20) admit an interesting interpretation in the context of 
the monodromy theory of the Fuchsian system of linear ODEs canonically related to the Riemann- 
Hilbert problem (2.6)-(2.9). We explain this connection in some detail in the Appendix. The 
results presented in the Appendix, however, are not used in the main body of the paper. 

To extend the theorem to nonzero V(z) we will use another differential identity. Let us introduce 
a parametric family of weights and the corresponding orthogonal polynomials indexed by t € [0, 1]. 
Namely, let 

(3.21) f(z,t) = (l-t + te v ^)e- v ^f(z). 

Thus f(z,0) corresponds to f(z) with V = 0, whereas f(z, 1) gives the function (1.2) we are 
interested in. 
Note that 

( ~ 99 , df(z,t) _ f(z,t)-f(z,0) 

Set now 7 = t and replace the function f(z) and the orthogonal polynomials in (3.5) by f(z, t) and 
the polynomials orthogonal w.r.t. f(z,t). Then the integral in the r.h.s. of (3.5) can be written as 
follows (we assume D n ^ for all n): 



(3.23) 



d_ 

di 



-f 



j7 / / .sSn{z 1 ,t) 2 I -1 .^ d 4>n(z,t) 
<j>n{Z,t) j 4>n(Z ,t)- 



2tt 



1 



Therefore, we obtain 



t (z,t) 
2n 



dz 

b„(z-\t) 



dz 



zf(z,t)d9 



dz 



-<f) n (z 1 ,t) 



l>n(z,t)A f( Z ,t)-f(z,0) 



_2n 1 f 
~T + 2^tJ c 



<f>n(z,t)- 



dz 

n(z'\t) 

dz 



- cj) n (z 1 ,t)- 



(z,t) 



dz 



dO 



zf(z, 0)dz. 



(3.24) ^]nD n (f(z,t)) = 2n[^ + \, 



'' dt 



2irtj c 



4>n(z,t) 



d<f> n (z 1 ,t) 

dz 



i{z~\t)- 



i>n{z,t) 
dz 



zf(z, 0)dz. 




Figure 1. Contour for the S-Riemann-Hilbert problem (m = 2). 



To write this identity in terms of the solution to the RHP (2.6) - (2.7), note first that using the 
recurrence relation (see, e.g., Lemma 2.2 of [15]) 

we have 

(3.25) Y 21 (z,t) = -Xn-lZ^n-liz- 1 ^) = - X nZ n Mz-\t) + M0,t)4> n ( Z ,t). 

Now using the orthogonality relations (2.1) and the formulae (3.2) and (3.3), we obtain from (3.24) 

Proposition 3.3. Let f(z,t) be given by (3.21) and D n (f(z,t)) ^ for all n. Let <j>k(z,t), <j>k(z,t), 
k = 0, 1, ... , be the corresponding orthogonal polynomials. Then for any n = 1, 2, . . . , 

<io^ °, n (ti + \\ 1 f -nf v i ^Y 21 (z,t) ^ dY u (z,t) \ df{z,t) A 

(3.26) -lnD n (f(z,t)) = — ^(z.t)—^— - y^^tJ—^—j ^^dz, 

where the integration is over the unit circle. 



4. ASYMPTOTICS FOR THE RlEMANN-HlLBERT PROBLEM 

The RHP of Section 2 was solved in [15]. In this section we list the results from [15] we need 
below for the proof of Theorem 1.1. We always assume (for the rest of the paper) that f(z) is given 
by (1.2) and that ctj ± f3j ^ —1, —2, ... for all j = 0, 1, . . . , m. In this section 4 we also assume for 
simplicity that zq = 1 is a singularity. However, the results trivially extend to the case ao = (3q = 0. 
In this section 4, we further assume that V{z) is analytic in a neighborhood of the unit circle. 

First, set 



(4.1) T(z) = Y{z) 



z~ na:i , \z\>l 
I, \z\ < 1. 



Now split the contour as shown in Figure 1. Set 



(4.2) 



S(z) 



T(z) 
T(z) 
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for z outside the lenses, 
for \z\ > 1 and inside the lenses, 
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1 

/(*) ' 
1 



1 ^— n ^ 







1 „n j 



for \z\ < 1 and inside the lenses. 



Here f{z) is the analytic continuation of f(z) off the unit circle into the inside of the lenses as 
discussed following (3.13). 

The function S(z) satisfies the following Riemann Hilbert problem: 

(a) S(z) is analytic for z € C \ E, where £ = U^ (E j U E^ U £"). 

(b) The boundary values of S(z) are related by the jump condition 



(4.3) 



S + (z) = S-(z) 



1 

f{z)- l z^ n 1 



zeu- o (£,u£;0, 



where the minus sign in the exponent is on Ej, and plus on 

/ n f(*\\ 

(4.4) S + (z) = S^(z 



/(*) 

(c) S{z) = 1 + 0(1/ z) as z -> oo, 

(d) As z — >■ Zj, j = 0, . . . , m, z <G C \ C outside the lenses, 



(4.5) 



(4.6) 



S(z) 



/o(i) o(i) + o(|z-^| 2a o 

\0{\) 0(l) + 0(\z- Zj \^) 



if «j ^ 0, and 



^.^(l) OQn\z- Zj \) 
* {Z) ~\0(1) 0(]n\z-Zj\) 



if a.j = 0, /3j ^ 0. The behavior of S(z) for z — > Zj in other sectors is obtained from these 
expressions by application of the appropriate jump conditions. 

We now present formulae for the parametrices which solve the model Riemann Hilbert problems 
outside the neighborhoods U Zj of the points Zj, and inside those neighborhoods, respectively. These 
parametrices match to the leading order in n on the boundaries of the neighborhoods U Zj , and this 
matching allows us to construct the asymptotic solution to the RHP for Y. 

The parametrix outside the U Zj 's is the following: 



(T>(z)°*, 



(4.7) 



N(z) 



V 
-1 0, 



\z\ > 1 
\z\ < 1 



V(z) a * 
where the Szego function 

(4.8) V{z) = exp — f l ^^ds, 

is analytic away from the unit circle, and we have 



zeC\uf =0 U z 



(4.9) 



V(z) = eH + (z) fj 



k=0 



z - z k 
z k e™ 



Olk+Pk 



\z\ < 1. 
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and 

m / y _ y, \ " a k+Pk 

(4.10) V{z) = b _ {z) -^l z _^\ , N> i, 

where Vb, b±(z) are denned in (1.8). Note that the branch of (z — z fc ) ±Qfc+/3fe in (4.9,4.10) is taken 
as discussed following equation (3.13) above. In (4.10) for any k, the cut of the root z~ ak+ @ k is the 
line 9 = 6k from z = to infinity, and 9^ < arg z < 2n + 9^. 

Inside each neighborhood U Zj the parametrix is given in terms of a confluent hypergeometric 
function. First, set 

(4.11) ( = nln-, 

z j 

where lnx > for x > 1, and has a cut on the negative half of the real axis. Under this trans- 
formation the neighborhood U Zj is mapped into a neighborhood of zero in the (-plane. Note that 
the transformation Q(z) is analytic, one-to-one, and it takes an arc of the unit circle to an interval 
of the imaginary axis. Let us now choose the exact form of the cuts £ in U Zj so that their images 
under the mapping C,{z) are straight lines (Figure 2). 

We add one more jump contour to £ in U Zj which is the pre-image of the real line T3 and T-j in 
the (-plane. This is needed below because of the non-analyticity of the function \z — Zj\ a i . Note 
that we can construct two different analytic continuations of this function off the unit circle to the 
pre-images of the upper and lower half (-plane, respectively. Namely, let 

(4.12) h aj {z) = \z-Zj\ a *, z = e w 

with the branches chosen as in (3.13). As remarked above, (3.13) is continued analytically to 
neighborhoods of the arcs < 9 < 9j, and 9j < 9 < 2tt. In U Zj , we extend these neighborhoods to 
the pre-images of the lower and upper half (-plane (intersected with ((U Zj )), respectively. The cut 
of h a . is along the contours T3 and T-j in the (-plane. 

For z — s- Zj, ( = n(z — Zj)/zj + 0((z — Zj) 2 ). We have < arg( < 2ir, which follows from the 
choice of arg(z — Zj) in (3.13). 
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From now on we will provide the formulae for the parametrix only in the region / (see [15] for 
the complete results). Set 

vu) m / z \ Pk/2 

(4.13) F j (z) = e— II - Y[K k {z)g h {z) 1/2 h aj {z) e ~ l ™ j i C € /, z€t^, j + 0. 

Note that this function is related to f(z) as follows: 

(4.14) Fj(z) 2 = f(z)e- 2 ^g^( Z ) ( € /. 

The functions gp k (z) are defined in (1.4). The formulae for Fq(z) are the same, but with gp {z) 
replaced with 



( 4 - 15 ) 5A,W = <^« n zeU Z0 . 



e -inf3 0j argz >0 
e ivr/3 0) argz<27r' 



We then have the following expression for the parametrix Pj(z) in the region z(I) of U Zj : 
(4.16) P*.(z) = Eiz^jiOFjizr^z^ 2 , ( e /. 



Here 



(4.17) £(z) = ^(z)^-^ ( ,),--/2 ^ q eftr(/ ,° +aaj) ) 



and 



_ / + ft, 1 + 2«„ C)e-( 2 ft'+^)e-C/2 

(4.18) vj/. ( g_ i _ c _ Qj ^ (1 _ q . + ^ 1 _ 2q . ? c)e , 7r(/3 ._ 3aj . )e - C / 2 r(i+a,+y 



< a ^(l + ftj - pj, 1 + 2aj, e - in C)e in(J,i+ai) e </2 T{ - +aj ~^ y 



ha,)„C/2l_ _ 

r(oy+&) 



where ip(a,b,x) is the confluent hypergeometric function of the second kind, and T(x) is Euler's 
T-function. Recall our assumption that ay ± fij ^ —1, —2, . . . 

The matching condition for the parametrices P Zj and N is the following for any k = 1,2,...: 

(4.19) P Zj (z)N' 1 (z) = I + A 1 (z) + A 2 (z) + --- + A k (z) + A ( [l 1 , z€dU Zj . 
Every A p (z), A p (z), p= 1,2,..., z 6 9U Z . is of the form 

(4.20) aj a3 0(n- p )af, aj = n^z~ n/2 . 

In particular, explicitly, on the part of dU z . whose £-image is in /, 
(4.21) 



which extends to a meromorphic function in a neighborhood of U Zj with a simple pole at z = Zj. 

The error term Aj^ in (4.19) is uniform in z on dU Zj . 
At the point Zj we have 

(4.22) Fj(z) = rije- 3i ™ j/2 z~ aj u a j (1 + O(u)), u = z- Zj , ( € /, 
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i out 




Figure 3. Contour T for the R and R Riemann-Hilbert problems (m = 2). 



where 
(4.23) 

and 
(4.24) 

(4.25) fij = [ e 



f 3-l 



,k=0 



k=j+l 



k+3 



Z k 



Afe/2 



(c^)) 2 = ^ {ai ~ Wi)n ~ 2Pj ( l + °(«)). « = C e J, 



&_(*;) 



Zj Zk 



The sums from to —1 for j = and from m + 1 to m for j = m are set to zero. 
4.1. R-RHP. Let 

\S{z)N~\z), zG[/oo\r, C/ oo = C\Uf =0 ^., 



(4.26) 



{ 1 \S(z)P-Mz), z€U Zi \T, j = 0,...,m. 



It is easy to verify that this function has jumps only on dU Zj , and the parts of Sj, S ■ lying outside 

the neighborhoods ?7 Z . (we denote these parts without the end-points X out , S out ). The full contour 
r is shown in Figure 3. Away from T, as a standard argument shows, R(z) is analytic. Moreover, 
we have: R(z) = I + 0(1/ z) as z — > oo. 
The jumps of i?(z) are as follows: 



(4.27) 
(4.28) 



R + (z) = R-{z)N(z)(^ f{z) \ z _ n f)N{z)-\ zG^\ 
= R-(z)N(z)( f ^_ lzn jW*)" 1 , z 



R+ z 



"out 
3 ' 



(4.29) R+{z) = R-(z)P Zj (z)N(z)~ 1 , z G \ { intersection points}, 

j = 0, . . . , m. 
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The jump matrix on £ out , S" out can be estimated uniformly in ctj, f3j as / + 0(exp(— en)), 
where e is a positive constant. The jump matrices on dU Zj admit a uniform expansion (4.19) in 

inverse powers of n conjugated by nP^^z- nCT3 / 2 5 and (4.20) is of order n 2max j l^ftl-p. To obtain 
the standard solution of the i?-RHP in terms of a Neuman series (see, e.g., [12]) we must have 
n 2ma Xj IStfil-i = (i^ that is fop. e (_i/2, 1/2) for all j = 0, 1, . . . ,m. However, it is possible to 
obtain the solution in any half-closed or open interval of length 1, i.e. for |||/3||| < 1, as follows. 
Let |||/?| 1 1 < 1 and consider the transformation 



(4.30) R(z) = n^ 3 R(z)n^ a3 z G C \ T, 
where 

(4.31) u = -(min^ + maxS^) 

2 j j 

which "shifts" all 3?/3j (in the conjugation terms of (4.20) for the jump matrix in (4.29)) into 
the interval (-1/2, 1/2). Note that u = 3?/3 io if only one $t/3 jo + 0, and u = if all = 0. 

Now in the RHP for R(z), the condition at infinity and the uniform exponential estimate / + 
0(exp(— en)) (with different e) of the jump matrices on S out , £ out is preserved, while the jump 
matrices on dU Zj have the form: 

(4.32) / + n wa3 Ax(z)n- wa3 + • • • + n^ 3 A k (z)n^ a3 + n^Al^zJn-^, z G dU Zj , 
where the order of each n^ 3 A p (z)n^ CT3 , n^ 3 A^ r) (z)n~ WCT3 , p = 1, 2, . . . , z G Uf =0 dU Zj is 

This implies that the standard analysis can be applied to the i?-RHP problem in the range G 
(q — 1/2, q + 1/2), j = 0, 1, . . . , m, for any q G R, and we obtain the asymptotic expansion 

k 

(4.33) R(z) = I + Y J Rp(z) + R$i(z), P = l,2... 

P =i 

uniformly for all z and for (3j in bounded sets of the strip q — 1/2 < 9?/3j < q + 1/2, j = 0, 1, . . . m, 
i.e. \\\/3\\\ < 1, provided aj ± /3j are outside neighborhoods of the points aj ± f3j = —1, —2, . . . (cf. 
(4.21)). 

The functions Rj(z) are computed recursively. We will need explicit expressions only for the 
first two. The first one is found from the conditions that R\{z) is analytic outside dU = U™ =0 dU Zj , 
Ri(z) — >■ as z — > oo, and 

(4.34) Ri,+ (z) = Ri,-{z) + n^ 3 Ai(z)n^ CT3 , z G <9f/. 
The solution is easily found. First denote 

(4.35) R k (z) = n'^ 3 R k (z)n^ 3 , R%\z) = n'"* 3 R ( k r) (z)n^ 3 , k = 1,2,..., 
and write for R: 

1 f A\(x)dx 



(4.36) Ri(z) = ±- [ 
2m J d 



du x 

m A 



A k (T \ i i m 77 

E^o^-A l( ,), *e^, ;=0,l,...,n,' 5C/ - U ^o^- 
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where the contours in the integral are traversed in the negative direction, and A k are the coefficients 
in the Laurent expansion of A±(z): 

(4.37) A 1 (z) = -^- + B k + 0(z-z k ), z -+ z k , fc = 0,l,...,m. 

z z k 

The coefficients are easy to compute using (4.19) and (4.25): 

The function i? 2 is now found from the conditions that R.2(z) — > as z — > oo, is analytic outside 
dU, and 



(4.39) R2,+ (z) = R2-(z) + R 1 -{z)n u<T3 Ai(z)n- u ' T3 + ra wa3 A 2 (z)n^ a3 , z G dU. 

The solution to this RHP is 

dx 
x — z 



(4.40) R 2 ( z ) = — f ( R 1 ^n^A^rT^ 3 + n WG ' A A 2 (x)n' ulC7 A — 

2vrz Jqjj V ' / x — 

At the fc'th step we have the RHP for R k (z) with the same analyticity condition and the condition 
at infinity, and the following jump: 

k 

(4.41) R kt+ (z) = R k ,-(z) + J2Rk- P A z ) n "' T3A p( z ) n ~ U " T3 , z e 9U, 

P =i 

where -Ro(z) = ^- 

We will now discuss the way in which the general R k (z) depends on n. In particular, we will 
discuss its order in n. First note that 

(4.43) bj = n^z~ n/ \ 5' ~ £V((A~A)-i) . 

Here the notation ^4 ~ B means A 1 = B', where X' is X in which each matrix element and 
each term in the sums is multiplied by a suitable constant independent of n. Starting with these 
expressions, and noting from (4.41) that 

k „ _ i_ 

(4.44) i? fc (z) - ^i? fc _ p ,_(z)n^A p (z)n-- CT3 - i^-i,-^)™^ 3 M*)™"^ 3 + -12 fc _i _(*), 

p=i n 
we obtain by induction: 

-2 



(4-45) i? 2p+1 (z) - ^tt ELo (^ 2 )* ( E ] 6 2 ) , 

(4.46) £ 2p+2 (z) ~ ^ ELo ^^2) - P = 0, 1, ■ ■ ■ • 



In particular, 

JJ 7 2 1 

'J "3 - / v "J - 1 ^ n 

(4.48) 0(5') = 0(5), 5 = max jifc n 2 ^"^)- 1 ) = ^(lllflll- 1 ), p = 0,l,..., 



(4.47) ^ p+1 (,) = ^o( E 1 fe2 ^V), ^ 2 (z) = ^0( 1 + ^ 2 
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as n — > oo. Here O(A) represent 2x2 matrices with elements of the corresponding order. 
Finally, note that the error term in (4.33) is 

(4.49) R^\z) = 0(\R k (z)\ + \R k+ i(z)\). 

In particular, as is clear from the above, if there is only one nonzero (3j , we obtain the expansion 
purely in inverse integer powers of n valid in fact for all (3j € C uniformly in bounded sets of the 
complex plane. 

It is clear from the construction and the properties of the asymptotic series of the confluent 
hypergeometric function that the error terms (z) are uniform for (3j in bounded subsets of the 
strip q — 1/2 < 3tf3j < q + 1/2, j = 0, 1, ... m, for aj in bounded sets of the half-plane Kay > —1/2, 
and for aj ± j3j away from neighborhoods of the negative integers. Moreover, the series (4.33) is 
differentiable in aj, j3j. 



5. ASYMPTOTICS FOR DIFFERENTIAL IDENTITIES AND INTEGRATION. PROOF OF THEOREM 1.1 



5.1. Pure Fisher-Hartwig singularities. The case ||/3||| < 1/2. First, we will prove the 
theorem for V{z) = and |||/3||| = max^fc \$t/3j — < 1/2. The proof is based on the differential 
identity (3.20). First, we show that (3.20) has the following asymptotic form. 

Proposition 5.1. Let (a®, f3o, . . . , a m , f3 m ) be in a compact subset, denote it A, belonging to the 
subset || |/3| || < 1, aj ± [3j ^ —1, —2, ... of the parameter space V = {(ao, (3q, . . . , a m , (3 m ) : aj, f3j G 
C, 5ftaj > —1/2} and including the point aj = j3j = 0, j = 0,1, ... ,m. Let f3j = if aj = 0, 
j = 0, 1, . . . , m, 5 = n^H^III -1 ). Then for n ->■ oo, and v = 0, 1, . . . , m, 



(5.1) 



d 



da u 
+ K + A,) 



\uD n {f{z)) = 2a v 



d T(l + a v + p v ) 

■ In — ; ; h m n 



da u ' r(l + 2a„) 



+ K - P„) 



d r(i + a„ - a,) 

■ In — : ; h m n 



da v r(l + 2a„) 



-E 

and 
(5.2) 



(a, + ft) In J— 51 + (aj-pj) In ■ 



2m^2(aj + Pj) + 0(n~ 1 \nn) + 0{5n\nn) 

3=0 







]nD n (f(z)) = -2(3 u 
d 



dPv 

+ (a„ + A, 



(ay + /3j) In ■ 



lnT(l + a„ + £„) - Inn 



(a,-/3,)ln^— 



+ (a, - /3„) 



9^ 



In r(l + a, -£„) + Inn 



v-l 



2iri ^2(aj + /?,■) + C^n" 1 In n) + 0{5n In n) . 



Remark 5.2. The error term 0(<5nlnn) = o(l) uniformly in A if |||/3||| < 1/2 — e, e > 0. In fact, 
the estimate for the error term can be considerably improved: see next section. 

Remark 5.3. The case (3j = if aj = is all we need below. After the proof of Proposition 5.1, we 
will integrate the identity (5.1) to obtain the asymptotics of D n for all aj / 0, (3j = 0. We then 
integrate the identity (5.2) and obtain the asymptotics of D n for all aj ^ 0, (3j ^ 0. This gives 
the general result for V = 0, \\\(3\\\ < 1/2, since we can set any aj = using the uniformity of the 
asymptotic expansion in the ay's. 
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Proof. Assume that for all j, (3j = if a 3 - = 0, and ^0, fe = 1,2,.... Then we can rewrite 

(3.20) in the form: 



r, dxn 
(5.3) ^i nAi(/(z) ) = _ 2 ^_ 



' + I 1 - n ( 2 n) (^4 n+ %) - ^r 1 ( r ) (^)y 2 ( r 1) (^)} 

3=0 

25S (?i ( 2 n) (^)^ 2 ( r 1) (^) -^^i ( r ) (^)? 2 ( 2 n+1) (^)) , 



3=0 

where 7 = ay or 7 = /3j. We now estimate the right-hand side of this identity as n — > 00. The 
asymptotics of Xn were found in ([15], Theorem 1.8). We need these asymptotics here in the case 
V = 0: 



^ m 

(5.4) x'-i = !--!>*-$ 



fc=0 



z ± Y n 2{h-^-i) v _i T ( l + <*j + ffj)r(i + Qfc - Pk) 

v k T( aj - Pj)T(a k + p k ) 
+ (5 2 ) + 0(5/n), 8 = n 2 ^~ 1 \ 



n — > 00, 



where 
(5.5) 



Vj = exp < 



j — 1 m 
-ITT [ ^«p- E 



a. 



n 



lz- - z l 2/3p 
^pl 



V P=0 p=j+l 

The asymptotics of y^(zj) were also found in [15] (Eq. (7.11)-(7.21)). Namely, 



(5.6) Y^{z ] ) = (I + rf ) )Lf\ L 



( n ) _ / M 2 \lijT}- 1 n OL i faz? M 2 2^jrjjn a i & 



'3 



where r 3 = pJ[\zj), the parameters rjj, fij are given by (4.23), (4.25), and 

^ r(i+a 3 -ft) r(2a ? ) \ 

M 



rCl+aj+^j) r(2aj) 
. r(l+2ey) n^j-Pj) . 



Note that the matrix has the structure 



(5.7) 



where L depends on n only via the oscillatory terms z™. 

Let us now obtain the asymptotics for the following combination appearing in (5.3): 

(5-8) 
d_ 



Y£\z 3 )Y^\z 3 ) = !-((! + r<g)L<g + rg>4?) (rr i} 4 2 +1) + (1 + rg^)^) 



c?7 



_ fr(n)r(n+l) n , Jn+lK , j (n) r (n+1) (n+l)\ 
- ^11 ^22 I 1 + r 22 J + ^11 ^12 r 21 J 

+ (4i } 42 +1) (i + + 4?^ST 1) r£ +1) 

We omit the lower index j of r and L for simplicity of notation. 



( + 11 ^ 11 + £7" 11 



r 12 ^^21 + ^r 12 



5 



97 

■(n) 



(n) 
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Now the explicit formula for L and the estimates for R(z) imply that 
(5.9) L^Lg +1) =0(l), = 0(1), 



(5.10) 
(5.11) 

where as before 
Therefore, 



T (n) T (n+l) (n+l) 
-^11 n 12 '21 



= 0(n5), 

^lnLH = 0(ln 



r (n) r (n+l) (n) 
^21 -^22 '12 

9 



0(n5), 



„(n+l) (n) 
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12 



0(<5), 



n J' ^r = 0(r) Inn, 
n 2(|||/3|||-i)_ 



(5.12) ^Y^\z 3 )y£ +1 \z 3 ) = L<S>LW £ in L W + Q ( — ) + In „ ) 



c^7 " 



c?7 



Inn 

n 



2aj z,- c?7 



Inn i . , 

+ C(dninn). 

n ' 



Similarly, we obtain 

(5.13) Y^\ Zj ) — ) 
and furthermore, 

(5.14) Y^\ Zj )Y^ +1 \z 3 ) = 0{5n) + 0(1), 



ln4i +1) + O ( — ) + 0(<5nlnn), 



Inn 

n 



FW(^)F 2 r ±; (^) = 0(5n) + 0(l). 



;(n+l) 



Note that because of the special structure of (5.7), the quantity n aj does not appear in any of the 
products (5.12)-(5.14). Substituting (5.12)-(5.14) into (5.3) and using the asymptotics (5.4), we 
obtain 



(5.15) ±-\*D n {f{z)) = ± 



3=0 



3=0 



+ ( — ) +0(5nlnn). 



n 



Let us calculate the logarithmic derivatives appearing in (5.15). From (5.6), (4.25), and (4.23) 
it is easy to obtain for the derivatives w.r.t. a u , v = 0, 1, . . . , m: 



(5.16) 
(5.17) 
(5.18) 
(5.19) 



-d_\ n L {n) - 9 In E0±^±M | i nn . 

da v mlj u,ll - da v m r(l+2a,) + mn , 



dav 3, 11 



(n) _ _ j n gj-z,. 



" InL 



+ 2vri, j < v; 
(n+i) 



9q 



r(i+2aj 

75^ hi ^li = 

OQi; J, 11 



a ln £fez^l +lnn; 



_ ln £2_i^ 



r(i+2aj) 



(n+l) 



3 ^ v. 



da v 3,21 2 „e' 

Similarly, we obtain for the derivatives w.r.t. j3 u : 
(5.20) J- In L % = J- In r(l + a v + - In n; 



( 5 - 21 ) 
(5.22) 

(5.23) 



In 



3/3. 



27ri, j < i/; 



J-lnL (n) 



In- 



J > ^ 



^7 ln L [T2i 1) 



InL 



d/3 v 
(n+l) 
J,21 



^-lnr(l + a„ -p v ) +lnn; 



J / v. 
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Combining these results with (5.15) we obtain (5.1), (5.2) on condition that Dk(f) ^ 0, k = 1, 2, ... , 
or equivalently (see Section 2), (ao, /?o, . . . , a m , (3 m ) £ A\f2fc . This condition can be replaced simply 
by («0) A)> • • • ) fi m ) G A in the following way. Let 

(5.24) /3 = «i = • • • = a m = /3 m = 0. 

Let f2fc (a;o) be the subset of fifc with a.,-, /3j fixed by (5.24). Since D^(f) = Dk(f(ao;z)) is an 
analytic function of «o and -Dfc(l) ^ 0, the set rifc ( a o) is finite. Let us rewrite the identity (5.1) with 
v = and assuming (5.24) in the form H'(ao) = 0, where H(ao) = D n (f(oto; z)) exp(— r(n, s)ds) 
and where r(n, ao) is the r.h.s. of (3.20) with 7 = ao and assuming (5.24). Since the expression 
(5.1) for r(n, ao) holds uniformly and is continuous for ao £ A provided n is larger than some /co(A), 
and D n (f(ao; z)) and its derivative are continuous, the function H {ao) is continuously differentiable 
for all n > ko(A). Hence, H'(ao) = for all ao € A and n > fco(A). Taking into account that 
H(0) = D n (l) 7^ 0, we conclude that D n (f(ao',z)) is nonzero, and that the identity (5.1) under 
(5.24) is, in fact, true for all ao £ A if n is sufficiently large (larger than ko(A)). Now fix ao £ A 
and assume the condition a± = ■ ■ ■ = a m = (3 m = 0. A similar argument as above then gives that 
D n (f(ao, /3o; z)) is nonzero and the identity (5.2) with v = is true for all /3o £ A if n is sufficiently 
large. Continuing this way, we complete the proof of Proposition 5.1 by induction. 

Remark 5.4. A similar argument applies to the asymptotic form of the differential identity (3.26) 
we need in Section 5.3 below. We omit the discussion. 

□ 



We will now complete the proof of Theorem 1.1 in the case |||/3||| < 1/2, V(z) = by integrating 
the identities of Proposition 5.1. In this case we denote 

D n {f{z)) = D n (a ,...,a m ;P ,...,p m ). 

First, set m = and /?o = 0. Then (5.1) becomes 

(5.25) A In D„ M = 2a. (l + ln„ + + <*.- In n) + 0(„„ln „). 

Integrating both sides over ao from to some ao and using the fact that D n (0) = 1, we obtain 
(5.26, D. M - „«|i±££ , 

where G(x) is Barnes G-function. To perform the integration we used the identity 



dxT(l + 2x)J G(l + 2z)' 

which easily follows from the standard formula (see, e.g. [28]): 

(5.28) I lnr(x + l)cfa; = -ln27r- ^ + ^ + z\nT(z + 1) - \nG(z + 1). 
Jo 2 2 

Now set m = 1, ao fixed. Set (3q = fi\ = 0. Relation (5.1) for v = 1 is then 

(5.29) -^-lnA t (a ,ai) = 2Q 1 ( 1 + Inn + -f- ^ + ^\ ) + 



ctai \ dao r(l + 2aij 

ao In zq + ao(ln z\ + i7r) — 2qq ln(zo — zi) + 27rzao + 0(n~ x Inn) + 0(5nlnn). 
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Integrating this over a\ from to some fixed a\ along a path lying in A (see Proposition 5.1) and 
using (5.27), we obtain 

(5.30) M "o.ai ) = a 2 lnn + 2i n G!(l + a 1 ) -lnG(l + 2ai) + a aiM^ie 37ri ) 
D n (a ,0) 

— 2ao«i ln(zo — z\) + 0(n _1 Inn) + 0(Sn Inn). 

Substituting here (5.26), we obtain 

(z - Z\) 



(5.31, D „ (ao , ai ,=n«3^jn|i^ 



zo^ie 3 ""* 



(1 + 0(n _1 In n) + 0(5n In n)) 



n Q o+ Q ? TT ^il±^l!| Zo _ Zl \- aoai (l + O(n^lnn) + 0(<5nlnn)) 
££G(l + 2ay) 

(to write the last equation we recall (3.13), the way the branch of (z — Zj) aj was fixed there, and 
the fact that argzo < argzi). 

Continuing this way, we finally obtain by induction for any fixed m, j3j = 0, j = 0, 1, . . . , m, the 
asymptotic expression 
(5.32) 

D n (a , ...,a rn ) = nP?=°°S + ^ JJ \ z . _ Zfc |-«^(l + O^ 1 Inn) + 0(<5nlnn)). 

3=0 ^ a i> 0<j<k<m 

We now add in the /3-singularities. We will make use of one more identity, which follows from 
(5.28): 



(5.33) / 
J o 



a + lnlYl + a + x) + (a — x) — lnlYl + a — x) — 2x I dx 

dx dx 



_ G(l + a + /3)G(1 + a - /3) 



G(l + a) 2 

First, we obtain the result for the case when — 1/4 < < 1/4. This implies that the order of the 
error term 0(5nlnn) remains o(l) as we integrate over /3j's starting at zero. As before, we always 
assume integration along a path in A. Setting v = 0, (3j = 0, j = 1, . . . , m in (5.2), and integrating 
this identity over /3q from zero to a fixed /3o, we obtain using (5.33): 

^ 5 34 ^ ln J B (tto,...,a m ;ft) = lnn + ln G(l + a + A))G(1 + a - /3 ) 



D„(a 0) ... ,a m ;0) " G(l + a ) 2 

+ /3 ^2 a j ln — h 0(n _1 ln n) + 0(<5n ln n) . 

Substituting here (5.32), we obtain 

(5.35, «,,(«„, ■ ■ ■ , «n>;A) = ^ «C^y ^-W U 3£S>1 

x ]j |«j - 2fcr a, ' a * IJ l^ - ) (1 + 0(^ 1 ln n) + 0(<5n Inn)). 

0<j<k<m j=l ^ ^ ' 

Next, set v = 1, f3j = 0, j = 2, . . . n in (5.2) and integrate over (3\. We obtain then the determinant 
D n (ao, . . . ,a m ; (3q, Pi). Continuing this procedure, we finally obtain by induction at step m the 
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asymptotics (1.10) for the case —1/4 < $tfij < 1/4 with V = and the error term 0(n _1 lnn) + 
0{5n\an) = o(l). 

Consider now the general case \\\/3\\\ < 1/2. We can choose q such that all $t/3j £ (q— 1/4, (7+1/4). 
Divide (0,g) into subintervals of length less than 1/2. Apply the above integration procedure to 
move all f3j from zero to the line where the real part of all f3j is the right end of the first subinterval. 
Since the length of subintervals is less than 1/2, the error term 0(<51nn) remains o(l). Recall that 
during the integration we avoid any points where ctj + j3j or ctj — (3j is a negative integer. Next, 
move the fij's to the right end of the second subinterval, and so on, until the point $tfij = q. From 
that point move the /3/s as needed. We thus obtain Theorem 1.1 with V = and \\\(3\\\ < 1/2. 

5.2. Pure Fisher-Hartwig singularities. Extension to \\\/3\\\ < 1. We now show that in fact 
the error term in (1.10) remains o(l) for the full range < 1. First, recall the definition of R 

and u in (4.30) and (4.31). We have \\\(3\\ \ = 2max i (5?/3 j - co) and 

-\<Wj-»<\ 

for all singular points Zj. We denote pj = Zj if $1/3 j — uj > 0, and m + the number of such points. 
Furthermore, denote qj = Zj if $1/3 j — oj < 0, and m_ the number of such points. Finally, let r,j = Zj 
if ). v \, co 0. 

Separating the main contributions in n (see (4.21)), we write the jump matrix for R on dU Zj (cf. 
(4.32)) in the form 

(5.36) I + n^A^n-^ + ■■■ = ! + A 1 (z)+D(z) + 0(n- 1 - p ), z£dU Zj , 
where 

(5.37) P= 1-|||)9||| 
and 

(5.38) D(z) = n^A^z)^^ - A\, z € dt/*., 

(5.39) Ai(z) = (n^ 3 Ai(z)n-^) 2 ia_ = -^^-<7_, cr_ - ' ° ° 



z-pj 

-2 



1 



M*) = - T V aj ~/^ (-^T-X p- n e-^-^^^n-^, z e dU p ., 



(5.40) A 1 (z) = (n^A 1 (z)n-^) 12 a + = ^-a + , <j + ~ <U ' 



z-qj 

2 







r ( a 3-Pj) \C PjF j( z )J n\n(z/ qj ) 

(5.41) Ai(z)=0, zedU rj . 

Note that -D(z) and Ai(z) are meromorphic functions in a neighborhood of U Zj with a simple 
pole at z = Zj. We have (see (4.24), (4.25), and recall that V = 0) 

(5.42) 6, = 6 W , hm W = -n^-^V ^i^f = 

(5.43) a,- = a< n) = lim a,-(z) = n -^-^)-i g "+i^ r(l + «j + ^) = Q(n _ p) ^ 

J J 7 r(ay - Pj) 
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Also note that 

(5.44) D(z) = 0(n- 1 ), zedU Zj . 

The main idea which will allow us to give the required estimate for the error term in (1.10) is 
the following. Write R in the form 

R(z) = Q(z)R(z), 

where R(z) is the solution to the RHP: 

(5.45) R(z) is analytic for z e C \ UjdU Zj 

(5.46) R(z)+ = R(z)-(I + Ai), z e UjdU gj , 

(5.47) R(z) = I + 0(l/z), z^oo. 

We solve this RHP below explicitly, however, note first that the solution exists and is unique 
and 

(5.48) R(z) = I + 0{n- p ) 

uniformly in z by standard arguments. We then have on dU Zj using (5.48), the jump condition 
(5.46), and the nilpotency of Ai: 

Q+ = R+R+ 1 = R-{I + Ai(z) + D(z) + OirT 1 ^))^ 1 = Q_(7 + D(z) + O^-P)). 

The jump matrix for Q on S out , S " out remains exponentially close to the identity. Therefore, 

D(s) 



-ds 



(5.49) Q(z) = I + Q 1 (z) + 0(n- 1 -<>), Q 1 (z) = ^- [ 

J 3 

In what follows, we will be interested in the matrix element Y2i(z) of our original RHP in a 
neighborhood of z = 0. In this neighborhood, we have using (5.49), 



(5.50) Y(z) = R(z)N(z) = n' UJa ' i Q(z)R(z)n UJrT3 N(z) = 



o r 
-i or 



From here, noting that V(0) = 1 as V(z) = 0, and using the estimates R(z) = I + 0(n p ), 
Qi(0) = 0(n _1 ), we obtain 

Xl-i = -Y 2 i(0) = (1 + Qi i22 (0) + O^- 1 -"))^^). 

Using the expression for Q\ from (5.49) and (5.38), we can write this equation in the form (recall 
that the contours dU Zj are oriented in the negative direction): 

(5.51) xl-i =(^~k " ® + j ^ 22 (°)- 

Eventually, we will use the product of these quantities over n to represent the determinant. Before 
doing that, we now solve the RHP for R and find i?22(0). 
Set 



(5.52) 



R{z) z eC\u jU z 

%)(/ + Ai(z)) z€UjU gj . 
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So defined, &(z) is obviously a meromorphic function with simple poles at the points Zj, which 
tends to I at infinity. Therefore, it can be written in the form 

m + $+ rn - $- 
(5-53) ^) = ^ + E^ + E^T 

3=1 P] 3=1 % 

for some constant matrices 

Moreover, the function R(z) (we use here the definition of Ai) in U Pk , k = 1, . . . , m+, 

(5.54) %) = $(z)(/ + A 1 )- 1 = 



'+E-^+E 



z-pj ' z-qj 

3=1 J 3=1 



z-pk 

is analytic in U Pk . Hence, the coefficients at negative powers of z — pt vanish. Equating the 
coefficient at (z — Pk)~ 2 to zero, we obtain 

(5.55) $+<r_ = 0, 

and therefore the matrix has the form 



(5.56) $ 



k 



9k 
fk o 



for some constants g^, fk- The vanishing of the coefficient at (z—pk) 1 gives the following condition 
on 3>^, where we used the equation 3>tcr_ = 0, 

(5.57) H~ h kY^ _ = ^ a -' fc = l,...,m+, 

J . =1 Pfc Qj 

with 6^ given by (5.42). Similarly, using the analyticity of &(z)(I + Ai(z))^ 1 in U qk , we obtain for 
all k = 1, . . . , m_ 



(5.58) <&- 



e k 

o %r 



(5.59) $ fc -a fc V^ = a fc cr+ 

^ Qk ~ Pj 

for some constants e&, Conditions (5.59), (5.57) are equations for the constants g k , fk, e&, hk- 
In view of (5.51), we are interested in 

(5.60) £ 22 (0) = 1>22(0) = 1-^^, 

3=1 qj 

where we used (5.53) to write the last equation. To calculate this quantity we first substitute <$>~j 
from (5.57) into (5.59). Then the 22 element of the resulting equations for k = 1, . . . , m_ can be 
written as follows: 

(5.61) (I-A)h = B, 

where h and B are m_-dimensional vectors with components and 

B k = Y] — k - J —, fc = l,...,m_, 
£i Ik ~ Pj 
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A = A^ is an m_ x m_ matrix with matrix elements 
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A k e = V akbj 
k/ JTi (Qk ~ Pj)(Pj -QtV 

and 7 is the m_ x m_ identity matrix. 

Define the m_ x m_ diagonal matrix A as follows 

A = diag{-<7i,-g 2 , • • • , -q m _}- 

Then (5.61) can be written in the form 

Tx = y, T = I - A^AA, x = A _1 /i, y = A- 1 B 

By Cramer's rule 

det(Ti"-y-T m _) 



dctr 



i=i 



detT 



where Tj are the columns of T and y is in the place of the fc'th column. We are interested in 
(5.62) 

First note that 

(5.63) detT = det(l - A _1 AA) = det(7 - A) = det(7 - A^). 

Second, a direct calculation shows that 



_ ^ _ det(7\ +y,T 2 + y,-- - ,T m _+y) 



(Ti + y,T 2 + !/,-•• ,T m _+y) = I- A', A' jk = 



a/ jb > 



\ (Qj ~Pe){Pe ~ QkY 



where 



a'j = ajqj \ 



b' £ = b e pe- 

Using the definitions (5.43,5.42) of aj, bj, we note that 

a r = af- l) + Oin- 1 -?), b'j = bf^ + Oin- 1 -?), 
and therefore A' = A^ 1 ^ + 0(n _1_p ). Thus we can rewrite (5.62) as 

(5.64) &,(„) - MO) - 1 - E f r [' + ° 

Note that 

det(J- A< n >) = l + 0(n- 2p ). 
Recalling now (5.51) and the representation of the Toeplitz determinant D n (f) as a product of 
Xu 2 -, we can ; f° r some sufficiently large no > 0, using (5.64), write 



n 



(5.65) D n (f) = D no (f) H Xk-i = 



fc=no+l 



Dno(f) n 

fc=nQ + l 



k i+ P 



det(7- AW) 
det(7- A^ 1 )) 



1 + 



C(n , a , . . . , am, A), ■ ■ ■ , /3m)n E ^=° ( ^ ^ 



1 + 



i 
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where all the error terms are uniform for ctj, f3j in compact sets, and C is a constant depending 
analytically on aj, (3j, and no only. Recall that in the derivation of this expression we assumed 
that \\\(3\\\ < 1 (and as usual !Raj > —1/2 for all j). Under this condition p > (see (5.37)), and 
the error term tends to zero as n — > oo. In the previous section, we obtained C explicitly for j3j 
satisfying < 1/2. Obviously, if all aj, (3j belong to fixed compact sets and the value of no 

is fixed, the constant C in (5.65) is bounded above in absolute value by a constant independent 
of any aj, (3j. By Vitali's theorem this fact implies that C can be analytically continued in (3j to 
the full domain |||/3||| < 1 off its values on the domain |||/3||| < 1/2. Thus C is given by the same 
expression also for |||/3||| < 1, and this concludes the proof of Theorem 1.1 for |||/3||| < 1, V(z) = 0, 
with the error term o(l) = 0(n'"^'" _1 ) in (1.10). 

5.3. Adding special analytic V(z). In this section and in the next one, we will add the multi- 
plicative factor e v<yZ \ where V is analytic in a neighborhood of the unit circle C, to a symbol with 
pure Fisher-Hartwig singularities and obtain the asymptotics of the corresponding determinant. 
Consider the deformation of the symbol f(z,t) given by (3.21) for t G [0, 1]. The analysis is based 
on integration of the differential identity (3.26) over t. In the present section we assume that V is 
such that 

(5.66) 1 - 1 + te v{z) + 0, t€[0,l], z£C, 

and 1 — t + te v ^ has no winding around C for all t £ [0, 1]. Then the Riemann-Hilbert problem 
(for the polynomials orthogonal) with f(z,t) has the same singularities as the problem with f(z) 
and is solved in the same way. In the following section we remove the condition (5.66). 

Let us rewrite the identity (3.26) in terms of the function S(z) = S(z, t) which is the solution 
of the Riemann-Hilbert problem posed on the deformed contour depicted in Figure 1. From the 
transformation Y — > S defined in (4.1), (4.2) we obtain 

(5.67) Y 11 = z n f- 1 S 12 , + + S 11 , + , Y 21 = z n f- 1 S 22 , + + S 21 , + , zeC = ^'. 

Substituting these expressions into (3.26) and taking into account that detS = 1 we arrive at the 
formula: 

(5.68) | In D n (f(z, t)) = n jf Z" 1 /^ + f Q [-/"V + (^2 >+ Sn,+ " S' 12 , + S 21 , + )] f^- 

+ J [ Z ~ H ( s 2i,+ s n,+ - S' lli+ S 2 i,+) + z n (S' 22+ S 12!+ - S' 12+ S 22}+ ) f~ 2 ] 

where we introduced the notation / = df/dt and /' = df /dz. Using the jump relation (4.4) 
satisfied by the function S(z, t) across the unit circle, we can rewrite equation (5.68) in a more 
symmetric way, 

f) r a 7 

(5.69) -lnD n (f(z,t)) = n f-^f +X (t), 

J c 

/i i i i i ■ dz 

^ [ s 22,+ s n,+ - S' 12t+ S 2 i, + + S' U _S 22 - - S 21 -S12-] f~ f— 

/■ dz 
[z~ n {S' 22 _S\ 2 ^ - S' 12 _S 22 -) + z n {S' 22+ Si 2 , + - S[ 2i+ S22,+)] f~ f^~- 

Analytically continuing the boundary values of S(z, t) from the " + " side of C to the " — " side 
of and from the " — " side of C to the " + " side of E, we can write the term X(t) in (5.69) in 
the form 

(5.70) X(t) = / (J. + «-)/-/^ + Jj-U + 



TOEPLITZ DETERMINANTS II 



27 



where 

(5.71) I = {S 22 S\ 2 - S' 12 S 22 )f \ J = S' 22 Sii - S' 12 S 2 i. 

Denote by S e (resp., £") the part of £ (resp., £") which lies inside WjL U Zj . Consider first 

Using (4.3) on £" and (5.71), we easily obtain that J_ = J + — z n I + and /_ = I + , and therefore, 

where C'J is the part of \Jp_QdU Zj lying inside the unit circle from the intersection of dU z . with 
the incoming £" to the intersection with the outgoing £" for each j. Note that J(z,t) has no 
nonintegrable singularity at 2^. Indeed, for z inside the smaller sector formed by £" at Zj, we can 
write 

(5.74) 5(z, t) = y(z, t) = Y(z, t) (I K{ \ ^ , 

where n(z,t) = c\(z,t){z — Zj) 2a \ if aj 7^ 0, and n{z,t) = c 2 (z, t)\n.{z — zj), if oij = 0, /3j 7^ 0, for 
some Ck(z, t) analytic near Zj, chosen so that Y{z) is analytic in a neighborhood of Zj. Writing J in 
terms of the matrix elements of Y, we see that the contributions of the (singular) derivative k'(z) 
cancel, and we obtain 

(5.75) J = {%Y n - Y 21 %)k + % 2 % x - %Y 2l . 

We now analyze (5.73) asymptotically. The asymptotic expression for S(z,t) inside the unit 
circle and outside UJL C/ Zj is given by (see (4.26), (4.7)) 

a g(z,t)* 3 ' 1 



n — > 00. 



(5.76) S(z,t) =R(z,t)e?W» J 

where g(z) is defined by the formula 

(5.77) V(z,t) = e 9 ^ l) . 
We have 

J = g' + R' n R22 ~ R12R21 =9 + O e (l/n), n -> 00, 

and we finally obtain 

(5-78) / {J_ + z-I_)r l f^.= I Jr 1 f^=[ 9'(z)r 1 f^ + £ (l/n), 
J-^/j 2m Jc'j 2m Jc'j 2m 

Similarly, using the jump condition for S and then the asymptotics for S outside the unit circle, 
we obtain 
(5.79) 

/ & <- U + = - I Jr'f^ = J a ^ + 0,(1/"). n =o, 

where C e is the part of \Jjt. dU Z:j lying outside the unit circle from the intersection of dU Zj with 
the incoming £ to the intersection with the outgoing £ for each j. 

Returning to the integrals (5.70), we now consider the part arising from the integration over 
£"\£" and £\£ £ . In this part the terms containing z ±n / T give a contribution which is exponentially 
small in n, while the integration of the terms with J can be replaced by the integration over £'\£^., 
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where T,' £ = £' n (U^L C/ 2j ), and over parts of the boundaries dU Zj . Thus, recalling also (5.78), 
(5.79), we have 
(5.80) 



n — ?■ oo, 



where d£/+ (resp., dU~) is the part of the boundary of U Zj inside (resp., outside) the unit circle 
oriented from the intersection with the incoming £' to the intersection with the outgoing £'. 

Note that using the same considerations as before, we can write in (5.80) J+ — J_ = g' + + cf_ + 
O e (l/n), and therefore 



(5.8i, x lt) - / »v/-v^ + / 9 '_/-v^ + Q,(i/«), 



n —f oo, 



where the closed anticlockwise oriented contours 

(5.82) £+ = (£' \ 4) u™ o at/+, s_ = (£' \ 4) u™ o 

(Note that one can deform £ + (resp., £_) to a circle around zero of radius 1 — e (resp., 1 + e).) 
By (3.21), 

(5.83) r 1 / = ~ 1 + e ^\ = I- In f 1 - t + te y . 
Furthermore, writing g in the form 

(5-84) g(z,t)=g Sz (z,t)+g FH (z), 
we have, by (5.77), (4.8), (4.9), (4.10), 

(5.85) g ^,t) = [ H1 - t + teV <">^, ,f» W _/».<««» + W<». 

Jc 2»r lEE=i(-«*+A)ln^, M>1 



Note that the solution of the Riemann-Hilbert problem is uniform in t G [0, 1]. From this fact 
and the explicit formulas (5.83), (5.84,5.85), we conclude, by an argument similar to the argument 
following (5.23) that the identity (5.69) holds for all t € [0, 1]. Now integrating equation (5.69) from 
t = to t = 1, we connect the Toeplitz determinant D n (f(z, 1)) with the Toeplitz determinant 
D n (f(z,0)) which represents the "pure" Fisher-Hartwig case and whose asymptotics we evaluated 
in the previous section. First, using (5.83) and changing the order of integration in the first term 
of (5.69), we obtain 

Furthermore, by (5.81), (5.84,5.85), 

(5.87) f dtX(t) = I Sz + I FH + £ (l/n), n -> oo, 

Jo 

where (cf. [17], Eq (86), (87)) 

(5.88) I s * = fdtl ((/% + (/*)'-) r 1 /^ = £ kV k V. k , 
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and 

(5.89) 7™ = (/«(,)); V W £ + / e (/»(-))'_ V W £ 

£ [ K + &) / mil + + A) / rm. - m\ * 

f^A Jn + z-z k 2m Jz_\z-Zk z J 2m 



fc=0 

Since 



(5.90) 9 Sz (z,l)= f 

Jc 



Sz . f V(s) ds \nb + (z) + V , \z\ < 1 



I c s — z2iri 1— ln6_(z), |z| > 1 ' 
we obtain 

(5.91) lnM*)=/ ln6_fe) = -/ 
which finally gives 

m 

(5.92) J™ = [-K + A) In &-(**) + (-a, + A) \nb + (z k )\ • 

fc=0 

Collecting (5.86), (5.87), (5.88), and (5.92), we obtain from (5.69) 



(5.93) hiD n (f(z,l)) -\nD n (f(z,Q)) = nV + Y, kV kV-k 

k=i 

m 

+ Yl H«fc + lnb -(*k) + {-oik + Pk) ^b + (z k )} + O e (l/n), n -> oo, 

fc=0 

which, in view of the result of the previous section, concludes the proof of Theorem 1.1 for analytic 
V(z) satisfying the condition (5.66). 

5.4. Extension to general analytic V(z). Now let V(z) be any function analytic in a neighbor- 
hood of the unite circle. Since zeros of the expression 1 — t + te v ( z \ t £ [0, 1], z € C, can only 
occur if 5fV(z) = ir(2k + 1), k G Z, there exists a positive integer q such that ^V(z) satisfies the 
condition (5.66) of the previous section, i.e., 

l-t + te« v(z) ^0, Vie [0,1], zeC, 

and this function has no winding around C for all t € [0, 1]. Let 

(5.94) fo(z) = f FH (z), ft(z) = e 1 < v Wf t - 1 (z), £ = l,...,q, 
where f FH (z) is the symbol for the "pure" Fisher-Hartwig case. Note that 

(5.95) f(z) = e v ^f FH (z) = f q (z). 
Consider fe(z), I = 1, . . . , q, and introduce the deformation, 

(5.96) f e (z, t) = (l - t + fe* v «) f^iz) = (l - t + te^) e^ V ^f FH (z). 
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All the considerations of the part of the previous section between equations (5.67) and (5.81) go 
through with f(z,t) replaced by fy(z,t) and we arrive at the formulae: 

8 f dz 

(5.97) j f In D n (f,(z, t)) = n J f^fe^ + X(t), 

X(t) = / e+ Afi-V^ + I + 0.<l/»), „ oo. 

where, as before, the closed anticlockwise oriented contours 

(5.98) £+ = (£' \ K) ur=o dU+, E_ = (£' \ T,' £ ) U™ dU~, 

and g(z) = ge(z) now corresponds to fa. 
The first term in (5.97) yields (cf. (5.86)) 

In order to evaluate X(t), we write g in the form 

(5.100) g(z,t) = g Sz (z,t)+g Sz (z)+g FH (z), 

where g FH (z) is the same as in (5.85), and 

,rm^ Sz, * f ln ( X -t + te* V(s) ) ds „ s £-1 f V(s) ds 

5.101 g Sz (z,t)= / g Sz (z) = / 

J c s - z 2m q J c s - z 2m 

Then we obtain 

(5.102) f dtX(t) = I Sz + I Sz + I FH + £ (l/n), n^oo, 
Jo 

where, up to the replacement V — >■ |V, the integrals I 5z and are the respective integrals from 
the previous section, i.e., 

j oo ^ m 

(5.103) / 5z = -^W_ fc , i- F ^ = -^[-K + /3 fe )ln6_(z fe ) + (-a fe + /3 fe )ln6 + (^)]. 

9 fc=l ^ fc=0 

The term J 52 in (5.102) is given by the equation 

(5.104) I Sz = I j c ({9 Sz {z))' + + i9 Sz {z))'_)v{z)^-.. 
Note that 

^ fc=i ^ fc=i 

Therefore, after a simple calculation we obtain 

(5.105) I s, = ±L_AY,kv k V- k . 

q k=l 
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Integrating (5.97) from t = to t = 1 and taking into account (5.99), (5.102), (5.103), and 

(5.105) , we obtain the following equation for the determinant D n (fy(z)): 

1 21 — 1 

(5.106) ln£> n (/*(*)) " lnf„(A-iW) = ~nV + —5- £ W_ fc 

q q k=i 

j m 

+ - V[-K + /3fc)ln6_(^ fe ) + (-a fc + /3 fc )hi6 + (z fc )] + O e (l/n), n -> 00. 
a 

^ fc=0 

This equation holds for any I = 1, . . . ,q. Summing up from I = 1 to t = q we again arrive at the 
formula 

00 

(5.107) In £>„(/(*)) - \nD n (f FH (z)) = nV + Y J ^kV-k 

k=i 

m 

+ [~(" fc + Pk) In6_(zfc) + (-a* + /3 fc ) In &+(**)] + O e (l/n), n -> 00, 

fc=0 

which conludes the proof of Theorem 1.1 in the case of V(z) analytic in a neighborhood of the unit 
circle. 

5.5. Extension to smooth V{z). If V(z) is just sufficiently smooth, in particular C°°, on the unit 
circle C so that (1.11) holds for s from zero up to and including some s > 0, we can approximate 
V(z) by trigonometric polynomials V^ n \z) = Yfk=-p{n) ^kZ h , z G C. First, consider the case when 
lll^lll = max^fc - $t.(3 k \ = 2maxj \$t/3j — uj\ < 1, where w is defined by (4.31). (The indices 
j, k = are omitted if «o = A) = 0.) We set 

(5.108) P=[n 1 -% v= III^HI +£i, 

where £1 > is chosen sufficiently small so that v < 1 (square brackets denote the integer part). 

First, we need to extend the RH analysis of the previous sections to symbols which depend on n, 
namely to the case when V in / is replaced by V^ n \ (We will denote such / by f(z, V^), and the 
original / by f(z, V).) We need to have a suitable estimate for the behavior of the error term in the 
asymptotics with n. For a fixed /, our analysis depended, in particular, on the fact that f(z)~ 1 z~ n 
is of order e~ £n , e' > 0, for z G S out (see Section 4.1), and similarly, f{z)~ 1 z n = 0{e~ en ) for 
z G S out . Here the contours S out , £ out are outside a fixed neighborhood of the unit circle (outside 
and inside C, respectively). If V is replaced by V^ n \ let us define the curve S outside U 1 JL U Zj by 



(5.109) 




7 >0, 



and S" outside WjL U Zj by 



(5.110) 




Inside all the sets C/ z . , the curves still go to Zj as discussed in Section 4. Let the radius of all U Zj 
be 2~flnp/p. We now fix the value of 7 as follows. Using the condition (1.11) we can write (here 
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and below c stands for various positive constants independent of n) 



(5.111) 



\z\ k 



\V^(z)\-\V \< £ \k s V k \ l -^-<c\ \ kSV ^ 
k=-p,k^0 1 1 



1/2 1/2 

{j^ (l_±3 7 lnp/p) ±2fc ^ 



< c 



K k=-p, k^O 

'A (l±3 7 lnfc/fc) ±2fc \ 1/2 



\k=i 



k 2s 



< C 



Y— 



\k=i 



fc2(s-3 7 ) 



1 + 



9 s n\ V 2 

ln 2 £: x n ' 



where z G £ out , z G dt/*. D {|z| > 1} (with "+" sign in "±"), and z G E out , z G n {|z| < 1} 
(with "— " sign). We now set 



(5.112) 
and then 



3 7 = a - (1 + e 2 )/2, e 2 >0, 



(5.113) |y (n) (z)|<c, |6 + (z,V (n) )| < c, |6_(z,F (n) )| < c, foralln 
uniformly on £ out , S" out , dl^.'s, and in fact in the whole annulus 1 - 3 7 ^ < \z\ < 1 + 3 7 ^. 

J P P 

It is easy to adapt the considerations of the previous sections to the present case, and we again 
obtain the expansion (4.19) for the jump matrix of R on dU Zj . Note that now \C( Z )\ = 0(n u \nn) 
and \z — zj\ = lnn/n 1 ~ u as n — > oo for z G dU Zj , and therefore using (4.19), (4.13), (3.13), (4.9) 
and the definition of v in (5.108), we obtain, in particular, 

(5.114) n^A 1 (z)n'^ = O (— ^— ) , z G Uf =0 8U Zj . 

\n £l Inn/ J J 

Furthermore, as follows from (5.109), (5.110), (5.113), and (4.27,4.28), the jump matrix on S out 
and S out is now the identity plus a function uniformly bounded in absolute value by 

(1— 1/\ 2maxj \St/3j\ / , \ T n 

^— j ii ±1 {i- v )^.\ <cexp{-|(l-^lnn}n 2 M m ^l^l, 

where the upper sign corresponds to S out , and the lower to S" out . 

The RH problem for R(z) (see Section 4.1) is therefore solvable, and we obtain R(z) as a series 
where the first term R\ is the same as before, and for the error term the same estimate holds for z 
outside a fixed neighborhood of the unit circle, e.g., for z large. 

This implies that Theorem 1.1 holds for f(z,V^). Note that it also holds for \ f(z,V^)\. 

We will now show that replacing with V in the symbol of the determinant D n (f(z, V^)) 
results, under a condition on s, in a small error only, so that Theorem 1.1 holds for D n (f(z, V)) as 
well. 

Using the Heine representation (2.10) for a Toeplitz determinant with (any) symbol f(z), the 
straightforward estimate 

uniformly for \z\ = 1, 



(5.116) 



b±(z,V^) = b±(z,V) 



l + O 



n 



1 



TOEPLITZ DETERMINANTS II 



33 



which follows from (1.11), and Theorem 1.1 for D n (\f(z, V^)\) and D n (f(z, F (n) )), we have if 
a(l -v) > 1, 



(5.117) D n (f(z, V)) — D n (f(z, V^] 



< 



(2vr) n n! 7 



2tt 



/•27T . 

. / jQ | e ^ _ e <*fc|2 JJ |/( e ^,y( n ))|d^ x ( 1 + c/nt 1 ^ 

l<j<fc<n J=0 

< ce 5RVon n E J "Lo((^) 2 +(^) 2 )( e c/™ (1 -" )s - 1 _ !) 



— i/)s 



< c 



Er = o((»«i) 2 +(Kft) 2 ). 



i 



< c 



A,(/(*,^ (n) ; 



n (l-i/)s-l 
-((l-z,) S -l-E"L ((^) 2 +(5R/3 3 ) 2 )) 



Therefore, 
(5.118) 

D n (/(^)) = D n (/(z,7<"')) | 1 + 



if 

(5.119) s > 



D n (f(z,V))-D n (f(z,V( n) )) 
D n (f(z,VW)) 



i + E7 = o((^) 2 + TO) 2 ) 

l-u 



= D n {f{zM n) ))(l + o{l)), 



Note that this condition is consistent with (5.112) and the requirement that 7 > 0. Using the 
expression for v in (5.108) and noting that E\ can be arbitrary close to zero, we replace (5.119) with 
(1.12). Under the condition (1.12) we then obtain the statement of the theorem for D n (f(z, V)). 



6. Appendix. The Toeplitz determinant D n as a tau-function 

In this section we construct a Fuchsian system of ODE's corresponding to the Riemann-Hilbert 
problem of Section 2 for V = 0. We show that the differential identities (3.20) for the Toeplitz 
determinant can be viewed as monodromy deformations of the tau-function associated with this 
Fuchsian system. 

Assume the pure Fisher-Hartwig case, V(z) = 0. Set 

m 

(6.1) $(z) = AY (n \z)A- 1 Y[{z- z k ) a ^z Xa \ 

k=0 

where 

k=0 k=0 

and the branches of all multi-valued functions are chosen as in Section 3. In terms of the function 
®(z), the Riemann-Hilbert problem (2.6)- (2.7) reads as follows: 

(a) $(z) is analytic for z G C \ (C U [0, 1] U {U^r,-}), where Tj is the ray 9 = 0j from zj 
to infinity. The unit circle C is oriented as before, counterclockwise, the segment [0, 1] is 
oriented from to 1, and the rays Tj are oriented towards infinity. 
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(b) The boundary values of 3>(z) are related by the jump conditions, 

,0 1 ' 



(6.2) $+(z) = $-(z) < 



z € C, 6j < argz < Oj+i, < j < m, 6 m +\ = 27r 



Z €Tj, 1 < j < m 
e - 27TiXa 3, z€[Q,l] 



where 

(6.3) Sj = exp < 



— ITT 



Oik 



k=0 



k=j+l 



k=0 



k=j+l 



(for j = m, the second and the fourth sums are absent) 
(c) $(z) has the following asymptotic behavior at infinity: 



3.4) 



*(*) = (i + o (£)) z^^+shfrt^+w 



as z — >■ oo and 0j < argz < < j < m, # m +i = 2ir (for j = m the last factor is 

omitted). 

(d) In the neighborhoods U Zj of the points Zj, j = 0, 1, . . . ,m, the function &(z) admits the 
following representations, which constitute a refinement of the estimates (2.8) and (2.9). 
• If aij / 0, then 



(6.5) 



*(z) = * j (z)(z-z j ) a i<»a 



3i 



where &j(z) is holomorphic at z = Zj (it is essentially the function Y(z) from Section 
3) and the matrix Cj is given by the formula, 



(6.6) 



with 



f l_ e 2T<(^+«j0 



(6.7) 



c i ~~ s j ' 



1 _ e 2 7r »(/3 j ~ Qj ) 



g4iricy l- e 2 "^J— 'j) 



47rja,- 5 

1— e J 



l Cj 
o iy 

zeU Zj , \z\<l 

z G f/ Zj , |z| > 1, argz < 9j 
z G C/ 2j ., |z| > 1, argz > 6>j 



(6.8) 



in the case j 7^ 0, and 



Co 



1 c 



5z > 



1 y j^Trm^ CJ Z <0 
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(6.9) 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



(6.14) 



with 



Co = so < 



' 1— e 27vi ^- l3 0+ a o) 
l_ e 2Tti(f3 -oi ) 



l-e 47ria 



o 47ria l-e 27r '^o-"o) 

^ e 47riaQ 



z G C4 , \z\ < 1 
2 e U Z0 , \z\ > 1, 9z < 
, z G C/ Zo , |z| > 1, > 



in the case j = 0. 
• If ay = and f3j ^ 0, then 

*(*) = *,-(z) (J ^ ln( x Z " Cj, d 3 = - Sj (l - e 2 ^), 

where &j(z) is again holomorphic at z = Zj and the matrix Cj this time is given by 
the formula 



1 Cj 

1 



with 



s,_i, z G J7 Z ., |z| < 1 



Cj = < 2 G C4,, |z| > 1, argz < 



dj, z G U Zj , \z\ > 1, argz > 0j 



in the case j ^ 0, and 
with 



1 c \ I /, 9z > 

1 J X 1 e 2^A CT3) q z < o 



Co = < 



do, 



2 G E/ 2o , |z| > 1, < , 



z G E/ 20 , |z| > 1, 9z > 



in the case j = 0. 

(e) In a small neighborhood of z = 0, the function $(2) admits a similar representation: 
(6.15) $(z) = ¥°\z)z Xa:i , 

where <&(°\z) is holomorphic at z = 0. 

We also note that all the matrices above have determinants equal to 1. 

A key feature of the $-RH problem is that all its jump matrices and the connection matrices Cj 
are piecewise constant in z. By standard arguments (see, e.g., [24] , [23]), based on the Liouville 
theorem, this fact implies that the function $(z) satisfies a linear matrix ODE of Fuchsian type: 



(6.16) 



d®(z) ., , „, . ^ A k B 
-A(z)$(z), A{z) = Y J — + - 



dz 



k=0 



z- z k z 
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with 

(6.17) B = A$( )(0) CT3 ($ (0) (0)) _1 

and 

a j ^ j {z j )a 3 ^~ 1 (z j ), if aj^O 



(6.18) A 



Moreover, since the jump matrices and the connection matrices Cj are all constant with respect to 
Zj, the function &(z) satisfies, in addition to (6.16), the equations 

(6.19) ^M = _^$( z) , j = l,...,m. 

OZj Z — Zj 

The compatibility condition of (6.16) and (6.19) yields the following nonlinear systems of ODEs on 
the matrix coefficients B and Ay. 

(6.20) eB = l A 11 B] t k+j 

(6 .21) dA L = _j^[A 11 M + [B 1 A j \^ 

dzj k=Q Zj — Zk Zj 

k^i 

In the context of the Fuchsian system (6.16), the jump matrices and the connection matrices 
Cj of the <1?-Riemann-Hilbert problem form the monodromy data of the system. The fact that 
these data do not depend on the parameters Zj means that the functions B = B(z\, ... ,z m and 
Aj = Aj(z±, . . . , z m ), j = 0, ...,m describe isomonodromy deformations of the system (6.16). 
Equations (6.20)— (6.21) are the classical Schlesinger equations. 

An important role in the modern theory of isomonodromy deformations is played by the notion 
of a r- function which was introduced by M. Jimbo, T. Miwa and K. Ueno in [24]. In the Fucshian 
case, the r-function is defined as follows. Let 

(6.22) Zj = J2^z=z k traceA(z)—^-^- 1 (z)dz k = ^ trace A k k ^ k) ^\z k )dz k . 

k=l k=l 
As shown in [24], the differential form Zj = Zj (Aq, . . . , A m , B;z±,..., z m ) is closed on the solutions 
of the Schlesinger system (6.20)-(6.21). The r-function is then defined as the exponential of the 
antiderivative of Zj, i.e., 

(6.23) ^ = lra ^%!$-( 2t ). 

It was already observed (see, e.g., [23] and [6]) that the r-functions evaluated on solutions of 
the Schlesinger equations generated by the Riemann-Hilbert problems associated with Toeplitz 
and Hankel determinants coincide, up to trivial factors, with the determinants themselves. In 
particular, for the Toeplitz determinant D n with the pure Fisher-Hartwig symbol, V(z) = 0, one 
can follow the calculations of [23] and obtain that 

(6.24) ^ = tr,ce A k ^pl^\z k ) + 2 £ - 2A^. 

dz k dz j^> Z 3~ Z k z k 
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Therefore, in the case of the Riemann-Hilbert problem (6.2)— (6.15), the relation between the asso- 
ciated Toeplitz determinant and the r-function is given by 



-2a j A 
Z 3 



(6.25) D n (f(z)) = D n ( Zl , z m \a,p) = r(zi, . . . , z m \a,p) H(z k - z 3 )~ 2a ^ J] . 

j<k j=0 

A direct substitution of formulae (6.5), (6.10), and (6.15) into equations (6.16) and (6.19) yields 

, d$ k (z k )~ i/ x trace AjAu 1 
6.26 trace A fe ^ $ k \z k ) = Y — + — trace BA k 

dz z k - Zj z k 

j^k 

= > trace A j —p^^- l (z j ) + traced — — ^ <£» (0) 0) . 

^ J dz k 3 y v dz k V V 

Hence equation (6.24) can be written as follows: 

Lemma 6.1. Let the Riemann-Hilbert problem for $ be solvable. For any k = 0, 1, . . . , m, 
(6.27) 

^ = Vtrace^^*-^) + trace B^ 0) (P»(0)) + 2 V - 2A^ 



3=0 i=° 



What is of interest here is that the differential identities (3.20), which play a very important 
role in the main text, can be written in a matrix form closely related to (6.27). For simplicity, we 
present only the case aj / 0. We have 

Lemma 6.2. Let the Riemann-Hilbert problem for be solvable. Let aj ^ ; j = 0, . . . , m. Then 
for any k = 0, 1, . . . , m, 

, dlnD n ^ A d^iizj)- , fl$(°)(0) /~ rnV 



i=o 



2 ay ln(zj — Zfe) — 2A ln(— Zfc) + ay In 2y — nlnz k , 



3=0 j=0 

&k 

OPk ~^ ORk OPk v j^o 

Remark 6.1. The significance of these equations is that they complement the isomonodromy defor- 
mation formula (6.27) by formulae which describe the monodromy deformations of the r-function 
(represented by the Toeplitz determinants D n ). (Equations (6.28) and (6.29) should be compared 
with the general constructions of the recent paper [7].) 

Proof. Although straightforward, it is rather tedious to derive (6.28) and (6.29) from the differential 
identities (3.20). There is, however, an alternative way to obtain (6.28) and (6.29) based on 
the direct analysis of the Riemann-Hilbert problem (2.6)-(2.7). First, we observe that the basic 
deformation formula for the Toeplitz determinant D n , i.e. equations (3.4), (3.5) can be written, 
using Eq. (2.4) in [15], in the form 

-on, 9 1 f _ n ( , dY 21 (z) dY n (z) v , A df(z) 
(6 - 30) = ^J C Z V n{2) — dz- Y2l{z) ) 
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Here, as in Section 3, 7 is either a k of j3k- Second, by 7-differentiating the Riemann-Hilbert problem 
(2.6)-(2.7) we easily obtain the following representation for the logarithmic derivative Y^ 1 (z) 
of its solution (cf. Lemma 2.1 of [7]): 

A /0 l\ v -i f ,,df(z') f ,,_ n dz> 



( , 31) XW = ^)y- W -^/y. M (S J)y. 



07 r — z 



1_ /■ ^-yn(z')F 21 (z') Y^fY) ^ dfW^-n dzf 



2mJc\ "W) ^11 CO^l (Z')J d 1 



z — z 



Now, from (6.1) and (6.16) we see that 

dY 11 (z] 



dz 



A u (z)Yn(z) + A u 2 A 12 (z)Y 21 (z) - c(z)Y u (z), 



= A 2 u A 21 (z)Y 11 (z) + A 22 {z)Y 21 {z) - c{z)Y 21 {z) 



where c(z) = YlT=o ~z=z^ + z - allows us to re- write (6.30) as follows: 
|- In D n = ±J z-nfatfA^Aniz) - A- 1 2 A 12 (z)Y 2 \(z) 

+ Y n (z)Y 21 (z) (A 22 (z) - A u {z)j) ^^dz 



A 2 Y^4 fv 2 (-\ df ^ Z ~ ndz A- 2 V^4 fv 2 (y\ df ^ Z ~ ndz 
An ^ A k>21 ^ Y U (z) — 2m{z _ - - A u ^ fc > 12 ^ y 2i (*) - 



A comparison with (6.31) yields the following /oca/ representation for the 7-derivative of lnl? n : 



(6.32) 



^ j / 1 /ft /ft /ft 

— lnD n = A^A fci21 X 12 (z fc ) + A n 2 ^A K12 X 2l {z k ) + A kA1 X u (z k ) + ^ A fe)22 X 22 (z fe ) 

fc=0 fc=0 fe=0 fc=0 

+ A 2 1J B 21 X 12 (0) + A n 2 5 12 X 21 (0) + Bnln(0) + B 22 X 22 (0). 



The last formula can be also written in the compact matrix form, 



rn 



(6.33) — ln£» n = trace A^AkAXizk) + trace A^BAA^O). 

^ fc=0 

By evaluating (6.33), with the help of the representation (6.5), one arrives at the formulae (6.28) 
and (6.29). □ 
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